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ABSTRACT 

Comparing the ejecta velocities at maximum brightness and narrow circumstcllar/interstellar Na D 
absorption line profiles of a sample of 23 Type la supernovae (SNe la), we determine that the prop- 
erties of SN la proge nitor systems and explosions are intimately connected. As demonstrated by 
Sternbe rg et al.l (|2011f) , half of all SNe la with detectable Na D absorption at the host-galaxy redshift 
in high-resolution spectroscopy have Na D line profiles with significant blueshifted absorption relative 
to the strongest absorption component, which indicates that a large fraction of SN la progenitor sys- 
tems have strong outflows. In this study, we find that SNe la with blueshifted circumstellar/interstellar 
absorption systematically have higher ejecta velocities and redder colors at maximum brightness rel- 
ative to the rest of the SN la population. This result is robust at a 98.9-99.8% confidence level, 
providing the first link between the progenitor systems and properties of the explosion. This finding 
is further evidence that the outflow scenario is the correct interpretation of the blueshifted Na D 
absorption, adding additional confirmation that some SNe la are produced from a single-degenerate 
progenitor channel. An additional implication is that either SN la progenitor systems have highly 
asymmetric outflows that are also aligned with the SN explosion or SNe la come from a variety of 
progenitor systems where SNe la from systems with strong outflows tend to have more kinetic energy 
per unit mass than those from systems with weak or no outflows. 
Subject headings: supernovae: general — distance scale — dust, extinction 



1. INTRODUCTION 

Type la supernovae (SNe la) appear to come from 
carbon-oxygen white d warfs (WDs) in a binary system 
(jHoyle fc Fowlerl [i960'): however, it is still unknown 
if they result from the mer ger of two WD s (the 
"double-degene rate" scenario; llben fc Tutukovi 119841 : 
[Webbmk|[l984) or a single WD accreting material from 
a non-de generate c ompani on ( the "singlc-dcgcncratc" 
scenari o; JWhelan fc Ibenl 119731: llben fc Tutukovi I1996T ). 
iHowelll ([201 ll ) and iWang fc Had (|2012l ) present re- 
cent reviews on the subject. With the exception 
of the recent l y discovered, nea r by SN la 20 life 



several different mo d els for SN la 
Nomoto et all Il984t iWooslev et all 



1991 



" Iffilie 
iPakn 



Hillcbrand t fc Niemeyerl 2000 : iGamezo et all 



explosions (e. 
19861: lKhokhlo\ 



2005[ iPakmor et al.l l2010t Ivan Kerkwiik et all 120101: 



Sim et all I2O1O0 and limited observational constraints 
(jBlondin et al.l 120111 ). This is particularly trou- 



dLi et all 120111: iNugent et all 120111: iBloom et all I201S. 
Brown et all 1201 lf Chomiuk et al. 2012; Margutt i etall 



bling since measurements of SNe la were used to 
discover the acc e lerating expansion of the Universe 
(jRiess et al.l 119981: iPerlmutter et al.l 11999ft and continue 
to be exqu isite tools for co nstraining cosmological 
parameters (iWood-Vasev et al.l 120071: iRiess et all 120071 : 



Freedman et aOT2009l: Ifficken et all| 2009b: K essler et all 



20091: lAmaiiu llah et al. 2010 MConlev etall 120111 : 



Suzuki et all 



2012). 



20121 ). there are only weak constraints on the na- 



ture of the binary companion for individual SNe la 
/e.g.. iMaoz fc Mannuccl l200l iFolev et al l [20101 [20ll 
lEdwards et all 120121: ISchaefer fc Pagnottal 120121) . and 

even then, it is difficult to completely exclude sce- 
narios for subsamples of SNe la. Similarly, there are 
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The recent observations of time variable, narrow (Av < 
100 kms -1 ; resulting from either circumstellar or in- 
terstellar material in the host galaxy) Na D absorp- 
tion in seve r al SNe l a (SNe 2006X 1999c l and 20071e; 
Patat et all l2007bl: iBlondin et all 120091: iSimon et all 



20091 ) are perhaps the best evidence that at least some 
SN la progenitor systems have strong outflows — the 
hallmark of certain classes of single-degenerate progen- 
itor systems. In the two cases with high-resolution 
spectroscopy, the variable component of the Na D ab- 
sorption profile was blueshifted relative to the strongest 
component of the profile. Using a large sample of 
SNe la with single-ep och high-resolution spectroscopy, 
ISternberg et al.l ([201 1L hereafter, Sll) determined that 
about half of all SNe la with Na D absorption have 
blueshifted Na D absorption relative to the position of 
the strongest absorption (12/22 SNe la). In the same 
sample, only a quarter of SNe la (5/22) have redshifted 
absorption, while the remaining quarter (5/22) have sin- 
gle or symmetric absorption profiles. If the absorption 
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were the result of only interstellar gas, we would expect 
the fra ction of blueshifted and redshifted profiles to be 
equal. IS 111 concluded that the difference in these num- 
bers was likely the result of a large fraction of SNe la 
(> 25% of SNe la with detected Na D absorption) com- 
ing from progenitor systems with strong outflows. 

The first SN la with ident ified variable Na D was 
SN 2006X (IPatat et alJl2007bD . This SN was also ex- 
ceptional in its hig h ejecta velocity and high reddening 
(jWang et al.ll2008| ). With only the detection of variable 
Na D in SN 2006X and the non-d etection of variable 
Na D in the low- velocity SN 2007af, iSimon et afl (|2007l ) 
hypothesized that variable Na D was exclusive to SNe la 
with high ejecta velocitie s. Since then, detect ions of vari- 
able Na D in SNe 1 999cl (|Blondin et al.ll2009D and 20071e 
(jSimon et a l. 2009) have been made — both SNe la with 
relatively high Si II A6355 velocities, «si n , near maxi- 
mum brightness. 

Additionally, ejecta velocity provides insight into the 
explosion. Directly, the ejecta velocity (up to viewing 
angle effects of asymmetric explosions) probes the ki- 
netic energy per unit mass. If one assumes that all 
SNe la have approx imately the same total mass (e.g., 
iMazzali et al.l l2007f) and relatively small asymmetries 
(e.g.. lLeonard et aLll2005l ). then ejecta velocity is a good 
tracer of the kinetic energy of the explosion. 

The strong correlation seen between the velocity gradi- 
ent of the Si II A6355 feature near maximum brightness 
of a SN la and the relat ive redshifting of f orbidden lines 
in its nebular spectrum (Maeda et al. 2010J) places strong 
constraints on explosion models. The early-time velocity 
measures the characteristics of the outer layers of the SN, 
while the late-time velocity measures the characteristics 
of the inner ejecta. The connection can be interpreted 
as the result of an asymmetric explosion, where different 
viewing an gles result in differ ent velocities at early and 
late times (jMaeda et al.|[20Toh . This asymmetric explo- 
sion interpretation can further explain the relation be- 
tween SN la ej ecta velocity (and velocity gradient) an d 
intrinsic color (Foley & Kascn 2011; Maeda et al. 2011). 

While ejecta velocity probes the characteristics of 
the SN explosion, the line profile of narrow Na D 
is influenced by the progenitor environment. There 
has been tenuous evidence that SNe la with variable 
Na D tend to have higher-velocity ejecta. To deter- 
mine if this trend is h appenstance, we use the large 
iFolev. Sanders, fc Kirshnerl (|201lL hereafter, FSK11) col- 
lection of maximum-brightness Si II A6355 velocities for 
SNe la, Wg; n , suppleme nted by additional data presented 
here to both expand the IS 111 sample and assign ejecta ve- 
locities to that expanded sample. 

The rest of this manuscript is structured as follows. 
In S ectio n [21 we present new data to expand the IFSK11I 
and IS 111 samples and then match the expanded samples. 
In Section [31 we examine the properties of the matched 
sample, finding that SNe la with blueshifted Na D line 
profiles tend to have higher ejecta velocities and redder 
colors than typical SNe la. We discuss this result and 
conclude in Section [4] 

2. DATA 

In this sect ion, we construct our sample. We begin 
with the IS 111 sample of SNe la and add to it additional 
SNe la with new high-resolution spectroscopy. We dis- 



cuss the IFSKIII sample of SNe la, and add SNe la with 
new photometry and low-resolution spectroscopy to the 
sampl. We exclude all SNe la with Amis(i?) > 1.5 mag 
since we cannot determine an accurate v^ n for these 
SNe. Our final sample, the F12 sample, is comprised of 
all spectroscopically normal SNe la with high-resolution 
spectroscopy and a measurement of Wgj n . 

2.1. High- Resolution Spectroscopy 

IS 111 present a sample of 35 SNe la with high-resolution 
spectroscopy primarily obtained with the Keck, Magel- 
lan, and VLT telescopes. For all spectra with any nar- 
row Na D absorption at the recession velocity of the host 
galaxy, a zero velocity for the system was assigned to the 
strongest Na D component. The SNe were then classified 
as having blueshifted, redshifted, single, or symmetric 
profiles. If, for example, the Na D profile had additional 
absorption offset primarily to the blue of the zero- velocity 
component, it was classified as Blueshifted. Depending 
on the relative strengths of circumstellar and interstellar 
absorption, it is possible in principle that some SNe la 
with strong-outflow progenitor systems could appear as 
Redshifted, Single, or (possibly) Symmetric. However, in 
a large sample of SNe with circumstellar material (CSM), 
the Blueshifted objects should outnu mber objects in the 
Redshifted sample. Of the 35 SNe la, HID find that 12, 5, 
and 5 SNe la have blueshifted, redshifted, and single or 
symmetric narrow Na D absorption profiles, respectively, 
while 13 have no detected Na D absorption. 

Using additional data, we modify the IS 111 sample. 
Specifically, we exclude SN 2008ge (which had no 
Na D absorption ) bec ause of its known peculiarities 
(|Folev et al.H2010D . IS 111 decided not to include some older 
SNe la in their sample because of a possible publication 
bias (i.e., SNe without detected or extraordinary Na D 
line profiles may not have ever been published). Since 
we are comparing velocities within our sample, this pub- 
lication bias should not affect us significantly. We there - 
fore add SN 1986G (IRichl 119871: iD'Odorico et al.lH98l . 
which has a blueshifte d Na D line profile, and SN 2001el 
(Sollerma n et al.ll2005h . which has a redshifted Na D line 
profile. Although SN 2000cx has high-resolution spectra 
(showing no Na D abs orpti on: IPatat et al.|[20~07al) . we do 
not add this SN to the IS 111 sample because of its peculiar 
nature (|Li et al.ll2001h . 

We also consider SN 1994D dHofc Filippenkd IJ995T : 
iKing et alJll995t IPatat et al.lll996f ). Its Na D line profile 
is complicated by its very low recession velocity. There 
are three groups of absorption lines at ~0, 200-250, and 
700 kms -1 , respectively. The feature at 700 kms -1 is 
clearly a single absorption line. If one interprets the fea- 
tures at ~200-250 kms -1 as being associated with the 
Milky Way, then SN 1994D would be classified as a Sin- 
gle object. However, if those features are associated with 
the SN or its host g alax y, then it would be classified as 
Blueshifted. In the IS11I sample, there were no features 
beyond ~150 kms -1 relative to the strongest absorp- 
tion feature, so a ^450-500 kms -1 offset would be an 
extreme outlier. Nonetheless, we do not claim to fully 
understand the characteristics of potential SN la progen- 
itor outflows. We consider the possibility of SN 1994D 
being a Blueshifted SN as low, and for our main analysis 
we classify it as Single. However, we also performed the 
analysis with SN 1994D classified as Blueshifted. The 



main results are not significantly affected by this choice. 

We add new high-resolution spectra of six SNe la 
obtained with the High- Resolution Spectrograph (HRS; 
Tull f 998) on the Hobby-Eberly Telescope (HET) in 200 f 
and 2002. The HRS observations used the GOO f/mrn 
grating centered at 6302 A, providing wavelength cover- 
age from 5300 to 7300 A at R = 30000, and were made 
through a 3" diameter fiber. 

The HET data were reduced in IRAF with the 
echelle package using standard procedures. Wave- 
length calibration was performed using spectra of a ThAr 
comparison lamp obtained immediately before or after 
the SN observations. The Na D emission lines in the 
HRS Hatfield lamps were masked out before flatfielding. 
Emission from the night sky was removed by subtracting 
the spectrum obtained through a sky fiber after scal- 
ing to match the sky emission line intensities in the SN 
spectrum (this scaling matches the relative throughputs 
of the science and sky fibers derived from internal flat- 
field frames when available). The normalized-flux spec- 
tra near the host-galaxy Na D doublet are presented in 
Figure [TJ 

We report the properties of the expanded high- 
resolution SN la sample in Tabled] 

2.2. Photometry 

IFSKllI presents 1630 v sj n measurements for 255 
SNe la ([Blondin eF al. 2012). Using a family of functions, 
they were able to estimate vsi n at maximum brightness, 
u Si id f° r ^Ne la with moderate decline rates and vg; n 
measured within about a week of maximum brightness. 
This large sample provides the reference sample for this 
work. Although there are potential differences between 
this sample and the underlying SN la sample, the size 
of the sample provides a reasonable baseline . To further 
increase the ove rlap between the expanded IS11I sample 
and the IFSKllI sample, here we provide new photome- 
tr y, spect roscopy, and/or derived values for seven SNe la. 

FSK11 presented several light curve-derived SN quan- 
tities (i max , -B max - Vmax, and Ami 5 (B)) that are crit- 
ical for measuring Wg s n . The light curves and original 
derive d quantities were part of the CfA3 dHicken et all 
2009al) and LOSS fc ancshalingam et alll2010D datasets 



Ganeshalingam et al.l (|2010f ) also presented light-curve 



param eters for SNF200805 14-002, which we use for this 
study. iFolev et al.1 (]2012t ) presented light-curve param- 
eters for SN 2009ig, which we also use here. For 
SN 2009nr, we use th e light-curve properties reported 
by iKhan et al.l (|2011l) . Light curves of SNe 2007on, 
2008C, 2008fp, 2 008hy , and 2008ia were presented in 
Stritzingc r et al.l ()2011D ; we derive light-curve parame- 
ters for these SNe by fitti ng the light curv es with the 
SNooPy Python package (jBurns et alll2011l ). We have 
corrected all measurements of B max — V max for Milky Way 
extinction. All light-curve parameters are presented in 
Table [2] 

Additionally, SNe 20091e and 2010A were followed by 
the FLWO 1.2 m telescope at Mt. Hopkins. Their light 
curves will be published as part of the CfA4 dataset 
(Hicken et al., submitted), but their light-curve param- 
eters are presented here. Again, the light curves were 
also fit with SNOoPy, and their light-curve parameters 
are listed in Table [2j 



2.3. Low- Resolution Spectroscopy 

In addition to the new light-curve parameters listed 
above, we present spectroscopic paramet ers here, specif- 
ically vsi ii and Wg; n . We use the IFolev et all (J2012D 
spectroscopic measurements for SN 2009ig. SN 2008C 
previously had sufficient light cur ves to determine a time 
of maximum brightness (jFSKllf ) but not a decline rate 
or B ma x — Knax- With additional photometric data, we 
were able to measure those quantities, but the measure- 
ment of tig; :l remains unchanged from that of IFSKllI 
For the remaining SNe, we discuss the data acquisition 
and reduction below. 

Additional low-resolution spectra were obtained with 
a variety of instruments on several telescopes. Specif- 
ically, we obtained d ata using the FAST spectrograph 
(jFabricant et al.lll998| ) on the FLWO 1.5 m telescope, the 
WFCCD spectrograph mounted on the Du Pont 2.5 m 
telescope, the IMACS spectrograph ()Dressler et al .1120061) 
mounted on the Magellan Baade 6.5 m telescope, and 
the LDSS-3 spectrograph (jAllington- Smith et al.l 119941 ) 
mounted on the Magellan Clay 6.5 m telescope. Specifics 
for each spectrum are presented in Table [2l 

Standard CCD processing and spectrum extraction of 
the low-resolution data were accomplished with IRAF. 
The data we re extracted using the optimal algorithm of 
Home (1986). Low-order polynomial fits to calibration- 
lamp spectra were used to establish the wavelength scale, 
and small adjustments derived from night-sky lines in 
the object frames were applied. The IMACS and LDSS3 
sp ectra were reduced w ith IRAF routines as described 
by lHamuv et al.l ([2006). For the remaining spectra, we 
employed our own IDL routines to flux calibrate the data 
and remove telluric lines using the well-exposed continua 
of spectrophotometric standards dWade fc Hornd fl988t 
IFolev et al.ll200llSilverman et al.ll2012D . 

All details of the low-resolution spectra are presented 
in Table H 

We determine «si i t and f§, T t fo r these S Ne using 
the same me thods of iBlondin et all (|2006D and |FSK11|. 
The IFSKllI method does not produce accurate Wg; n 
measurements for SNe la with Ami^(B) > 1.5 mag. 
As a result, we do not report Vg { u for SNe la with 
Amis(-B) > 1.5 mag. Our measurements for vgi n and 
u Si ii are re P or ted in Table [5J 

We present the full cross-matched sample of SNe la 
with both Wgj :I measurements and high-resolution spec- 
troscopy in Table [T] The SNe with both Ugj u measure- 
ments and high-resolution spectroscopy are designated 
as the 'F12' sample. 

3. RESULTS 

Using the F12 sample, we are able to compare dif- 
fere nt subsa mples to each other as well as to the F12 
and IFSKllI samples. There are currently only four 
SNe la in the Single/Symmetric subsample (when in- 
cluding SN 1994D). Although this subsample is small, 
it is large enough to perform basic statistical analyses. 
Nonetheless, we caution against over-interpretation of 
the results for the Single/Symmetric subsample, and sug- 
gest that additional analysis should be performed when 
the sample sizes increase. For the results presented here, 
we classify SN 1994D as Single; if SN 1994D is classified 
instead as Blueshifted, the significance of our results is 
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TABLE 1 

Properties of the Expanded High-Resolution SN Ia Sample 



SN 



UT Date 



Phasc a 



Ami 5 (B) 
(mag) 



iax Vm 

(mag) 



u Si II 
(10 3 kms" 1 ) 



Ambiguous 11 



Host 
Type 



Blueshitted 


1986G 


1986 May 9.4 b 


-1.9 (0.1) 


1.76 


0.01) 


0.88 


0.09) 






SO pec 


2002bo 


2002 Mar. 21.3 


-2.4 (0.7) 


1.24 


0.05) 


0.40 


0.04) 


-13.1 


0.2 


Sa 


2002ha 


2002 Nov. 2.6 C 


0.2 (0.7) 


1.35 


0.05) 


-0.09 


0.03) 


-11.3 


0.2 N 


Sab 


2002jg 


2002 Nov. 30. l d 


-1.4 (0.1) 


1.54 


0.04) 


0.61 


0.03) 






Pair 


2006X 


2006 Feb. 18 e 


-1.6 (0.8) 


1.18 


0.03) 


1.26 


0.05) 


-16.0 


0.2 


Sbc 


20071e 


2007 Oct. 20.4 


-5.4 (0.4) 


1.05 


0.04) 


0.29 


0.04) 


-13.3 


0.2 


Sc 


2008C 


2008 Jan. 17.4 


16.0 (0.3) 


1.07 


0.03) 


0.24 


0.03) 


-11.3 


0.2 


SO/a 


2008dt 


2008 Jul. 6.5 


6.5 (3.1) 


1.02 


0.14) 


0.55 


0.08) 


-14.3 


0.3 


SO/a 


2008fp 


2008 Sep. 17.3 


-4.1 (0.1) 


0.92 


0.01) 


0.49 


0.01) 


-10.7 


0.2 


SO pec 


2009ig 


2009 Oct. 16.3 


39.9 (0.1) 


0.89 


0.02) 


0.14 


0.04) 


-13.5 


0.2 


Sa 


Redshifted 


2001cp 


2001 Jun. 28.3 


0.1 (0.4) 


0.92 


0.04) 


-0.03 


0.03) 


-10.9 


0.2 


Sbc 


2001da 


2001 Jul. 22.4 


4.8 (1.3) 


1.25 


0.05) 


0.16 


0.03) 


-11.5 


0.2 


Sab 


2001el 


2001 Sep. 24.8 f 


-5.3 (0.1) 


1.15 


0.03) 


0.03 


0.03) 


-11.7 


0.2 N 


Scd 


2007af 


2007 Apr. 8.4 


24.2 (0.4) 


1.23 


0.05) 


0.00 


0.03) 


-11.0 


0.2 


Scd 


20091e 


2009 Nov. 26.0 


-4.8 (0.1) 


0.96 


0.02) 


0.12 


0.02) 


-12.1 


0.2 


Sbc 


2010A 


2010 Jan. 8.0 


-8.2 (0.1) 


0.88 


0.02) 


0.16 


0.02) 


-10.3 


0.2 


Sab 


Single/Symmetric 


2002cr 


2002 May 13.2 


-1.3 (0.5) 


1.23 


0.03) 


-0.05 


0.03) 


-10.1 


0.2 


Scd 


2006cm 


2006 Jun. 6.5 


7.0 (0.5) 


1.04 


0.07) 


1.07 


0.06) 


-11.5 


0.2 


Sb 


2007sr 


2008 Jan. 17.6 


33.3 (0.3) 


1.07 


0.07) 


0.09 


0.06) 


-12.9 


0.3, 


Sm pec 


No Absorption 


2006eu 


2006 Sep. 10.3 


2.3 (2.2) 


1.37 


0.14) 


0.49 


0.04) 


-11.4 


0.3 


E 


2007hj 


2007 Oct. 20.3 


43.8 (0.2) 


1.97 


0.06) 


0.12 


0.03) 






SO 


2007on 


2008 Jan. 17.3 


62.0 (0.1) 


1.64 


0.01) 


0.12 


0.01) 






E 


SNF20080514-002 


2008 Jun. 13.4 


17.0 (0.5) 


1.39 


0.04) 


-0.16 


0.05) 


-10.8 


0.2 


SO 


2008hv 


2008 Dec. 15.3 


-1.9 (0.1) 


1.33 


0.03) 


0.03 


0.01) 


-10.8 


0.2 


SO 


2008ia 


2008 Dec. 15.2 


2.0 (0.1) 


1.35 


0.01) 


0.00 


0.01) 


-11.5 


0.2 


SO 


2009nr 


2010 Jan. 9.3 


12.5 (0.3) 


0.93 


0.03) 


0.00 


0.03) 


-9.6 


0.2, 


Scd 



1994D 



1994 Apr. 14 1 



23.9 (0.4) 1.42 (0.01) -0.04 (0.04) -11.2 (0.2) 



SO' 



Note. — Uncertainties arc listed in parentheses. 
a Days since B maximum for high-resolution spectrum. 

Exposure-time weighted mean observation time from[D'Odorico ct al. (1989), which presented the highest S/N spec- 
trum. Observations span 1986 May 7.2 - 11.1. 
c Exposure-time weighted mean observation time. Observations span 2002 Oct. 30 — Nov. 6. 

Exposure-time weighted mean observation time. Observations span 2002 Nov. 28 — Dec. 2. 
e Date corresponding to the first epoch of spectroscopy from Patat ct al. (2007b) 

Exposure-time weighted mean observation time. Observations span 2001 Sep. 21.2 - 28.3. 
g Exposure-time weighted mean observation time. Observations span 2001 Oct. 16 — 17. 

See text for discussion of the classification of SN 1994D. 
1 Date for Ho £z Filippcnko (1995) observation, which presented the highest S/N spectrum. 
J The host galaxy shows an obvious dust lane in images. 



TABLE 2 
Log of Low- Resolution Optical Spectral Observations and Basic Parameters 



SN 


UT Date 


Instrument 


Time of 
Max. (MJD) 


Phase a 


Amis(B) 
(mag) 


'-'Si II 
(10 3 kms" 1 ) 


, Si II . 
(10 3 kms" 1 ) 


2007on 

2008C b 

SNF20080514-002 C 

2008fp 

2008hv 

2008ia 

2009ig d 

20091c 

2009nr c 

2010A 


2007 Nov. 17.2 

2008 Jan. 5.3 
2008 May 28.3 
2008 Sep. 21.4 
2008 Dec. 18.3 

2008 Dec. 15.3 

2009 Sep. 5.4 

2009 Nov. 24.2 

2010 Jan. 7.5 
2010 Jan. 16.1 


IMACS 

FAST 

FAST 

WFCCD 

LDSS-3 

LDSS-3 

LRS 

FAST 

FAST 

FAST 


54419.9 (0.1) 
54466.1 (0.3) 

54613.0 (0.5) 
54730.5 (0.1) 

54817.1 (0.1) 

54813.2 (0.1) 
55080.0 (0.1) 
55165.9 (0.1) 
55192.7 (0.3) 
55212.4 (0.1) 


1.3 (0.1) 

4.1 (0.3) 

1.3 (0.5) 

-0.2 (0.1) 

1.1 (0.1) 

2.1 (0.1) 

-0.4 (0.1) 

-6.6 (0.1) 

10.6 (0.3) 

-0.3 (0.1) 
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Fig. 1.— High-resolution HET spectra of SNe 2001cp, 2001da, 2002bo, 2002cr, 2002ha, and 2002jg. The spectra are presented in 
normalized flux units and only the region near Na D is shown. The data are presented on a velocity scale where the zero velocity 
corresponds to the strongest absorption. The red and blue spectra are the individual Na Dl and Na D2 components, respectively, while 
the black spectra represent the combination of the two components. We label each SN and its profile designation in each panel. 
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reduced slightly, but not significantly 

3.1. Phase 

The standard model to explain the variable Na D seen 
in SNe la requires ionizing photons from the time of ex- 
plosion (or immediately after) to io nize the Na atoms i n 
the circumstellar environment (e.g. JPatat et a l. 2007b). 
After sufficient time, the Na atoms would recombine, re- 
sulting in increased Na D absorption. If this model is cor- 
rect, one might worry that systems with outflows will on 
average have no preferred line profile at early times (when 
the Na is ionized), but will on average have blueshifted 
profiles at later times. 

To test this idea, we compared the different subsamples 
of the F12 sample. The F12 sample has a median phase 
of 1.0 days relative to B maximum. The Blueshifted, 
Redshifted, Single/Symmetric, and No Absorption sub- 
samples have median phases of —1.5, —2.3, 15.5, and 
12.5 days, respectively. There is no indication that SNe 
in the Blueshifted subsample tend to have later phases. 

There are three possible scenarios that allow for the 
Blueshifted subsample to predominantly come from pro- 
genitor systems with outflows and not have a phase bias. 
First, it is possible that the systems with variable Na D 
are well tuned such that the ionizing flux and the dis- 
tribution of the CSM are matched to show variability, 
while other systems with outflows simply do not ionize 
the CSM because of a lack of ionizing photons, the CSM 
being too far/dense, or both. Alternatively, it may be 
the case that the recombination timescale for most sys- 
tems is shorter than the median time between explosion 
and the time of the spectrum, or about two weeks (al- 
though this is unlikely since variable Na D is not seen 
until after maximum brightness even when there are ear- 
lier observations). Finally, progenitor systems may have 
Na at various radii from the SN, where the closest Na 
is ionized, but the more distant Na produces blueshifted 
absorption. Indeed, the two SNe la with high-resolution 
spectroscopy and variable Na D (SNe 2006X and 20071e) 
have blueshifted Na D absorption at all spectroscopic 
epochs. 

Regardless of why there is not an obvious phase bias in 
the F12 sample, it does not appear that the Blueshifted 
sample is missing a significant number of potential out- 
flow systems specifically because the Na in the CSM was 
fully ionized at the time of spectroscopy. 

3.2. Light- Curve Shape 

We would like to compare the decline rate distributions 
of the various samples. The full high-resolution sample, 
including all SNe la regardless of decline rate (i.e., the 
F12 sample combined with SNe la with high-resolution 
spectroscopy and Am\^{B) > 1.5 mag), has a weighted 
mean Am 15 (B) of (Ami 5 (B)) = 1.22 ± 0.01 mag. The 
Blueshifted, Redshifted, Single/Symmetric, and No Ab- 
sorption subsamples of the full high-resolution sample 
have weighted means of (Arn^i?)) = 1.21 ± 0.01, 
1.10 ± 0.02, 1.23 ± 0.02, and 1.45 ± 0.03 mag, respec- 
tively. The subsample of SNe la with no Na D absorp- 
tion have significantly faster light curves than the bulk 
of the sample. Six of 7 SNe la in this subsample come 
from SO or earlier host galaxies (Table [lj , which is not 
surprising since these host galaxies are expected to have 



significantly less gas and dust, while no SN hosted in 
an e lliptical galaxy is in any of the other subsamples 
(|S11[ ). SNe la in elliptical gal axies tend to have signifi- 
cantly faster light curves (e.g., Howell 2001; Hickcn et al. 
I2009bl) 

The iFSKllI method of determining Vg { u does not ap- 
ply to SNe la with Ami^(B) > 1.5 mag. As a re- 
sult, SNe la with Ami^(B) > 1.5 mag can not be in- 
cluded in the F12 sample, and are therefore not included 
in any further analysis. Since many of the SNe la in 
the No Absorption subsample of the full high-resolution 
sample have large decline rates, the F12 sample does 
not include many SNe la without Na D absorption. 
With this restriction on light-curve shape, the F12 sam- 
ple has a weighted mean Ami§(B) of (Ami§(B)) = 
1.14 ± 0.01 mag, and the Blueshifted, Redshifted, Sin- 
gle/Symmetric, and No Absorption subsamples have 
weighted means of (Ami 5 (B)) = 1.12 ±0.01, 1.10±0.02, 
1.23 ± 0.02, and 1.25 ± 0.02 mag, respectively. The No 
Absorption subsample is still skewed to faster declining 
SNe relative to the entire F12 sample. 

3.3. Maximum-Light Color 

We now compare the B max — V m&K pseudo-color (after 
correcting for M ilky Way reddening) for the various sam- 
ples. The iFSKllj sample has a median pseudo-color of 
0.01 mag. This is significantly bluer than the F12 sam- 
ple, which has a median pseudo-color of i? m ax — Knax = 
0.12 mag. A K-S test results in a p-value of only 0.048, 
indicating that the two samples may be drawn from the 
different parent population. 

We interpret the color difference as a sel ection effect. 
Since the SN sample from which the lFSKUl sample is de- 
rived is the sample of low-z SNe discovered in SN surveys, 
there is some amount of Malmquist bias associated with 
both the discovery and follow-up of the sample. This 
bias will naturally remove extinguished SNe la near the 
limit of SN searches from the sample. However, for high- 
resolution spectroscopy, SNe typically need to be much 
brighter than the magnitude limit of SN searches; as a 
result, the Malmquist bia s in the F12 sample should be 
significantly less than the IFSKllI sample. 

The median pseudo-colors for the Blueshifted, Red- 
shifted, Single/Symmetric, and No Absorption sub- 
samples are i? max — V max — 0.34, 0.08, 0.02, and 
0.00 mag, respectively. The values for the Blueshifted, 
Single/Symmetric, and No Absorption subsamples are 
fairly different from the median value for the full F12 
sample. However, the Single/Symmetric and No Ab- 
sorption subs amples have median pseudo-colors similar 
to that of the IFSKllI sample, perhaps further indicating 
that SNe la in the Single/Symmetric and No Absorption 
subsamples have minimal dust extinction. Performing a 
K-S test between the Blueshifted subsample and the F12 
sample excluding the Blueshifted subsample, we find a p- 
value of 0.021. A similar test for the Single/Symmetric 
and No Absorption subsamples results in p- values of 0.56 
and 0.16, respectively. Although the median values for 
all three subsamples suggest that there might be differ- 
ences in the underlying populations of the subsamples, 
only the difference between the Blueshifted subsample 
and the remainder of the F12 sample is statistically sig- 
nificant. 

The redder colors for the Blueshifted subsample are 
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Fig. 2.— CDF of u§. u for the IFSK11I sample excluding SNe 
in the F12 sample (black line, labeled 'FSK11'), the F12 sample 
(gray line, labeled 'F12'), and the Blueshifted subsample (blue line, 
labeled 'Blue'). 



consistent with there being some SNe la in the sample 
with CSM resulting in additional reddening that is com- 
bined with the underlying host-galaxy reddening distri- 
bution. 

3.4. Ejecta Velocity 

With the knowledge that Blueshifted and No Absorp- 
tion SNe la are redder and faster decliners than the rest 
of the F12 sample, respectively, while the other subsam- 
ples do not appear to be significantly different from each 
other, we now examine correlations between narrow Na D 
profile properties and Vg { n . Figure [H displays the cumu- 
lative distribution function (CDF) of v^ n for various 
sampl es. T he gray and black lines represent the F12 
and IFSKTll samples, respectively. We also display the 
Blueshifted subsample of the F12 sample with a blue 
line. 

The CDF of the F 12 sample has a very similar shape 
to that of the iFSKTll sa mple. Th e median Vg { n for both 
the F12 sample and the IFSKTll samples is —11,300 and 
—11,400 km s -1 , respectively. Performing a K-S test (ex- 
cluding any SNe in the F12 sample from the [FSKlj] sam- 
ple), we find a K-S p- value of 0.87. From t his test, there 
is no indication that the F12 and iFSKllI samples have 
different parent populations. 

However, the Blueshif ted subsample is qualitatively 
different from the lFSKlll and F12 samples. Interestingly, 
the five SNe la with the highest Vg ; n in the F12 sample 
are all in the Blueshifted subsample. Specifically, 63% of 
the Blueshifted subsample has a higher v^ n than any 
SN la in any of the other subsamples. Not surprisingly, 
the CDF of the Blueshifted subsa mple app ears to be sig- 
nificantly different from that of the FSK11 sample or the 
F12 sample, with the Blueshifted subsample tending to 
have higher v^ u than the typical SN la. The Blueshifted 
subsample has a median v^ n of —13,200 kms -1 , while 



the Redshifted, Single/Symmetric, and No Absorption 
subsamples have median Wg ; n of —11,300, —11,300, and 
— 10,800 kms -1 , respectively . A K- S test to compare 
the Blueshifted sample to the IFSKllI sample (again, ex- 
cluding any SNe la from the F12 sample) results in a 
p-value of 0.034, indicating that the two samples likely 
have different parent distributions. 

Although the F12 sample has a similar v^ { n distribu- 
tion as the FSK11 sample, there may be other important 
ways in which the parent populations of the two samples 
are different. To mitigate any potential bias, we per- 
formed K-S tests between the Blueshifted subsample and 
the F12 sample (excluding the Blueshifted subsamples). 
This test resulted in a p- value of 0.018. The Blueshifted 
subsample appears to be drawn from a different parent 
population than the complementary F12 sample. 

We tested the robustness of this result using a Monte 
Carlo simulation. For each iteration, we selected a sam- 
ple of 23 SNe (to match the number of SNe la for 
whic h we h ave high-resolution spectroscopy) from the 
full IFSKTTI sample of SNe la. From this sample, we then 
selected a subsample of "Blueshifted" SNe equal to the 
number of Blueshifted SNe la in the F12 sample. We 
performed a K-S test on these random samples similar 
to that described above. We also noted if the five highest 
or lowest velocity SNe were in the "Blueshifted" sample. 

We found that 1.1% of all random samples had a K-S 
p-value as small or smaller than that of the Blueshifted 
sample. This value is similar to the p- value that we found 
above (0.018), and indicates that the p- value is a reason- 
able measurement of the probability of the Blueshifted 
and F12 samples being randomly drawn from the same 
parent distribution. We also found that 0.28% of the ran- 
dom samples had the five highest or lowest velocity SNe 
in the "Blueshifted" subsample. It is important to test 
both extremes since either case would have been notable 
to us. This test of the extreme SNe is perhaps a more in- 
teresting statistic since it is highly plausible that some of 
the Blueshifted systems come from the same parent pop- 
ulation as the Redshifted systems, but are simply con- 
fused through random line-of-sight effects. Furthermore, 
only 0.21% of the random samples had both of the above 
characteristics. Although this last number is something 
of an a posteriori measurement, it does indicate that the 
number of extreme objects in the Blueshifted sample is 
providing additional information beyond the simple K-S 
statistic. 

Compared to the F12 sample, the Blueshifted subsam- 
ple appears to have a statistically significant difference in 
Ug; n . As seen in Section |3~U1 the SNe in the Blueshifted 
subsam ple also tend to have r edder observed B maK — V maK 
colors. iFolev fc Kasenl (|2011ft and IFSKllI showed that 
ejecta velocity and intrinsic color are highly correlated 
with SNe la with higher ejecta velocity tending to be in- 
trinsically redder. In Figure [31 we sho w v §, TT a s a func- 
tion of observed -B max — Kn ax for t he IFSKllI and F12 
samples. We also display the 1FSK 1 ll relation between in- 
trinsic -B m ax — Knax and Wg ; n . Several SNe la (including 
the majority of the B lueshifted subsample) fall far off to 
the red of the IFSKllI relation, indicating that their very 
red colors are dominated by dust reddening rather than 
the trend with intrinsic color. 

In the color-velocity parameter space, the Blueshifted 
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- V m ax pseudo-color for the 
IFSK11I (black crosses), Blueshifted (blue circles), Redshifted (red 
squares), Single/Symmetric (purple diamonds), and No Absorption 
subsamples (gold triangles). The solid line is the velocity-color 
relation between «g ; n and intrinsic B max — Vmax IjFSKlll l. The 
empty black circles indicate SNe la with detected variable Na D. 
The region define d by the dashed lines (corresponding to a 2-<r 
deviation from the FSK11 relation between intrinsic B max — V ma) . 
and i)g ; I: and iig ; n = —12,000 kms -1 ) has a large fraction of 
SNe la with blueshifted and/or variable Na D line profiles. 



subsample is even more distinct from the rest of the Ff 2 
sample. Performing a two-dimensional K-S test on the 
Blueshifted subsample, the resulting p-value of 0.023 is 
further indication of the Blueshifted subsample being a 
distinct subpopulation in this parameter space. 

In Figure [HJ we have also indicated the SNe la for 
which variable Na D has been detected. All three have 
high Wg ; n and i? max — Vmax- Two of those three have a 
blueshifted Na D profile, while the third does not have 
the high-resolution spectroscopy necessary to make that 
measurement. If one wants to choose SNe la which have 
a high probability of having variable Na D, one might 
want to choose reddened SNe la with high ejecta veloc- 
ities. We defined a region of this parameter space cor- 
responding to high velocity (wg ; n < —12,000 kms -1 ) 
and red (redder than the 2-cr deviation from the IFSK11I 
relation between intrinsic -B max — V„ lax and Ug; n , which 
roughly corresponds to 5 max — V max > 0.2 mag). This 
region of parameter space contains half of the Blueshifted 
subsample and all SNe la with observed Na D variability. 
We suggest that SNe la in this region are excellent can- 
didates for having variable Na D. Additionally, special 
attention should be paid to SNe that have Blueshifted 
Na D absorption in the first epoch. 

4. DISCUSSION & CONCLUSIONS 

For the first time, we have shown a direct connection 
between the progenitor system of a SN la and its ex- 
plosion properties. In contrast, previously, only indi- 
rect connections between progenitors and SN la prop- 
erties have been made by examining host galaxy proper- 



ties (e .g., lHamuv et aT1ll996l lHowellll200l1: iHicken et all 
l2009bf) . Specifically, we find that SNe la with blueshifted 
narrow Na D profiles tend to have higher velocity ejecta 
as probed by Si II A6355 at maximum brightness than 
those with no Na D absorption or those with redshifted, 
single, or s ymme tric profiles. This was achieved by ex- 
panding the lSUl sample of SNe la with additional SNe la 
with high-resolution spectroscopy and matchin g that ex- 
panded high-resolution sample to the iFSKllI sa mple of 
SNe la with measured Wg ; n , again, expanding the lFSKill 
sample slightly with new data presented here. The re- 
sulting sample is designated the 'F12' sample. 

This result is further evidence that a significant frac- 
tion of SN la progenitor systems have outflows (|S11[ ). and 
that the large fraction of SNe la with blueshifted Na D 
line profiles is not caused by some other effect. Similarly, 
this result supports the interpretati on that the variable 
Na D is the resul t of ionized CSM (jPatat et al.|l2007bfc 
Simon et a l. 2009) rath er than a geometric effect (jChugail 
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For most subsamples and the larger SN la population, 
there is no indication of different light-curve shape distri- 
butions. However, the No Absorption subsample tends 
to have faster declining light curves (as expected given 
the host-galaxy distribution of the subsample). There is 
no phase bias for the Blueshifted subsample, indicating 
that few systems with outflows are artificially placed in 
other subsamples because the circumstellar Na is com- 
pletely ionized at the time of spectroscopy. 

The F12 sample has on ave rage redder i? max — VJ nax 
pseudo-colors than the lFSKill sample. We suggest that 
this difference is the result of Ma lmquist bias reducing 
the number of reddened SNe in the lFSKill sample. Com- 
pared to the rest of the F12 sample, the Blueshifted sub- 
sample is slightly redder and the No Absorption sub- 
sample is slightly bluer. The bluer color for the No Ab- 
sorption subsample is consistent with a SN population 
that has minimal dust reddening. Similarly, the redder 
color for the Blueshifted subsample is consistent with a 
SN population with additional dust in the circumstellar 
environment. 

High-resolution spectra of SNe la in E/S0 galaxies have 
never shown Na D absorption (with the exception of 
SNe la hosted in peculiar SO galaxies, such SN 1986G, 
which was in Cen A) . Most single-degenerate progenitor 
scenarios should exist in old stellar populations, so the 
different host properties for the different subsamples in- 
dicate that the population of SN progenitors with strong 
outflows may be extremely rare in old stellar populations. 
Alternatively, SNe la from early and late-type galaxies 
may have, on average, very different progenitor systems. 
Future studies of SN hosts and locations may constrain 
the age of these progenitor systems, further constraining 
the exact progenitor systems. 

There are clear implications of the connection between 
blueshifted Na D absorption profiles and high ejecta ve- 
locities. Specifically, if all SNe la come from a single 
progenitor channel, then one requires an asymmetric pro- 
genitor system and explosion where the higher-velocity 
ejecta are aligned with the higher density regions of the 
CSM or explosions with higher kinetic energy per unit 
mass also have a denser and/or closer CSM. Alterna- 
tively, SNe la could come from a variety of progenitor 
channels where those that originate in progenitor sys- 



terns with strong outflows tend to have more kinetic en- 
ergy per unit mass than those with weak or no outflows. 
From several different methods, specific progenitor sys- 
tems which shoul d have strong outflows have been ruled 
out for SN 201 lfe ()Li et al.ll201ltlBloom et alJl20ll . and 
there is strong eviden ce that its progenitor system had 
a clean environment (iChomiuk et al.ll2012t iPatat et al.l 
1201 It [Margutt i et al.l I2012D . Since a significant fraction 
of SNc la do have strong progenitor outflows, it appears 
that there are at least two progenitor channels for SNe la. 
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